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Abstract 

We consider the features of multiparticle tree cross sections in scalar theories in the framework of 
a semiclassical approach. These cross sections at large multiplicities have exponential form, and the 
properties of the exponent in different regimes are discussed. 
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Considerable interest has been attracted in recent years to the issue of multiparticle production both in 
perturbative and non-perturbative regimes in weakly coupled scalar field theory (for a review, see 0]). 
This problem has been initiated by the qualitative observation ^) that in the ordinary theory with the 
action 

s= ld^ x (\(d»ipf-\^-±y) a) 

(hereafter we set the mass of the boson equal to unity) cross sections of the processes of creation of a 
large number of bosons by few initial ones exhibit factorial dependence on the multiplicity of the final 
state. The reason is that the number of tree graphs contributing to the amplitude to produce n particles 
ly-j grows as «!. At n ~ 1/A this factor is sufficient to compensate the suppression due to the small coupling 

^\ • constant, and the tree level multiparticle cross sections become large. 

Now there is a lot of perturbative results (see, for instance, Q|) which confirm the factorial growth of 
the tree level amplitudes, and, what is more important, exhibit the exponential behavior of multiparticle 
cross sections 

Qh / 1 \ 

dJ . a(E,n) cx exp ( — F(Xn, e) ) (2) 

where e — (E — n) /n, and E is the energy of initial state. Though practically all these results have been 
obtained in the perturbative regime Xn <C 1, £ <C 1, one expects that the exponential behavior survives 
in the regime A — > 0, with An, e being fixed. Moreover, there are strong indications [|| that the exponent 
jj^ ■ F{\n,e) is universal, i.e. independent of the few-particle initial state. 

The exponential form of <r(E, n) implies that for the calculation of the function F(Xn, e) there may 
exist a semiclassical method like the Landau technique || for calculating matrix elements in quantum 
mechanics. In fact, some approaches aiming to generalize the Landau method have been proposed recently 
ill- 

Let us concentrate on the technique proposed in [Q . It is based on singular solutions of some classical 
boundary value problem. As the result, the function F(Xn,e) obtained by perturbative calculations 
has been reproduced by means of this way. In the (d + l)-space-time it looks like 

F(Xn,e) = Xnln^-Xn+^(hi— + l)+^-^-Xne + BX 2 n 2 + 0(X 3 n 3 ) + 0(X 2 n 2 e) + 0(Xne 2 ) (3) 
16 2 V na / 12 

where B is some known numericalconstant. Moreover, certain new properties of F have been obtained 
in different regimes. The given method is very promising since it allows one to study F(Xn, e) by the 
perturbative expansion or, in any case, to put this issue in the computer. 



lr To appear in Proc. Xth Int. Workshop on high energy physics and quantum field theory, Zvenigorod, 1995. 
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In this paper we study the behavior of F(Xn, e) in the regime An <C 1, with e being fixed, i.e. in fact 
we will investigate only the energy dependence of the tree level cross sections. We present first few terms 
in the low energy expansion of F in powers of s (up to Xne 2 ). A lower bound on the function F at any 
energies is obtained by numerical calculations. This bound seems to be larger than the estimate found 
by Voloshin |10| from the analysis of the Feynman diagrams and coincides with the asymptotics at the 
high energies [g. We demonstrate also that in fact the function F lies above this bound and does not 
coincide with it. 

The paper is organized as follows. In Sect. 2 we present a calculation of the boundary value problem 
for the tree level cross sections suggested in [[llj leading to the result which coincides with one obtained 
from the more general approach M. In Sect. 3 we find the 0(e 2 )-correction. Sect. 4 is devoted to the proof 
that, at least in high energy limit, spherically symmetric solutions do not provide the true value of F. 
In Sect. 5 some numerical results are presented. In Sect. 6 we present our conclusions. Several technical 
details are collected in Appendix. 



2 General formalism 

Now with the aim to introduce our notations and to expose some points of the calculation we will briefly 
discuss the method suggested in |ll[] . We will examine the cross section of the decay of one virtual particle 
at rest with the energy E into n real particles in the model ([j]). Let us start from the coherent state 
representation (see [jl2| for details) for the matrix element (P\Sip\0). In general, it has the form 

(p\Scp\0) = lim f DipDip f ip(E, k) exp{iS + B f (/3* k , ip f )} 

Tj -,oo J 

where 

p(25,k)= J dtd d xif(t,x)e~* Et +^, ip k (t) = J ^L_^, x )e- ikx , <^(k) = ip k (T f ) 

B f {Pt,v } ) = -\ j 'd d mi3*_ k e 2 ^ -~ Jd d ku ipf (k) ipf (~k) + Jd d kV2^e^ T ff3* kl p f (~k), u = Vl + V 

To calculate this matrix element at the tree level one should extremize (iS + Bf) over if and iff. This 
extremization yields the classical field equation with the boundary conditions 

dlip + if + Xip 3 = 0, ip k {t^ +oo) = -^=e*"*, fk(t -» -oo) ~a fc e Mt (4) 

V 2cj 



The solution to this equation, ip c ([3* , x) has only positive frequency parts due to the boundary conditions, 
and hence the action S and the boundary term Bf are reduced to zero. Thus, the tree level matrix element 
is 

F{(3*) = ({p}\Sf\0) tICC = J dtd d xip c ((3*,t^)e^ Et +^ 
and the amplitude 1 — > n is given by 



{/3*}=0 



This result coincides with one obtained in jl3j from LSZ-reduction formula. 

According to the coherent state formalism the n-particle cross section can be easily found, 



^.trcc 



d£ f n a 7-1 o* j_ _ la* _M f :t\I± xi\ i __/\ 1 



D(3D(3*dxdx'if c ((3*,x)if* c (P, x') exp { -iE(t - t') + ik(x - x') - - / d d k(3*(3 - nln^ 



(5) 
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where the integral over £ corresponds to the n-th derivative due to Cauchy theorem. In order to estimate 
the integral (||) the saddle point technique can be applied. To perform it correctly, however, one should 
look more carefully at ip c and get rid of zero modes in ip c corresponding to time shifts. One writes 

I3* k = b* k e^ ip c {p\t) = ip(b*,t + r,) 

where rj is the collective coordinate, and b k are new integration variables which obey a constraint to be 
spacified below. (Of course, there are zero modes corresponding to the space shifts, but we do not write 
them explicitly since taking them into account leads only to spacial momentum conservation). So, Eq.(||) 
becomes 



tree 
^l-m 



j^j J DbDb*dr 1 dri'(p(b*,E,k)ip*{b,E,k)exp^-iE(r 1 -r/) - ~J d d kb*be iu{ v' ~v) _ nln^j 



Now, if one imposes a constraint that ip does not contain exponents and, so, may be considered as a 
pre-exponential factor, then this integral can be calculated by making use of the saddle point technique. 
The exact form of this constraint, which may be viewed as a constraint on b k , will be determined in what 
follows. Therefore, the tree cross section is 

a(E, n) trco oc e^-" (6) 

where W cx tr is the saddle point value of the functional 

W = ET-n9- [ t ^kb%b k e u,T - e (7) 



with respect to T — i(jf — 77), 9 = ln£, bk, and b k under an additional constraint which is yet to be 
defined. 

In general, one expects that these classical solutions, being considered in Euclidean space-time, have 
singularities on some surface r = r(x). The leading behavior of the integral 

tp(E,k) = J d d+1 xip(x)e- tEt+lkx 

is determined by the nearest to the origin singularity of ip{x) in the region — r = Im t < 0. If this 
singularity is located at r = r*, then 



ip(E,k) 



oc e 



,-Et 



So, the constraint that ip may be considered as a pre-exponential factor requires r* = 0. Since we are 
working in the center-of-mass frame, the result should depend on |k|, but not on its direction. So, in 
x-space we require 0(3) spherical symmetry and finally the constraint is that ip{x) has the singularity at 
r = 0, x = 0. 

Therefore, the problem of finding the cross section at any value of E and n can be formulated in 
Euclidean space-time and consists of the following steps, 

• One selects from all solutions </?(t, x) to the Euclidean field equation, (9(3)-symmetrical ones which 
are singular at r = x = and have the following asymptotics at r — > 00, 

Ad u* 

(^rF^ (T ' X)e (8) 
From them one finds Fourier components b k and then determines W from Eq.(j7[). 

• One should extremize W over all bk, b* k (that is over different singular surfaces), T, and 9. The tree 
level cross section is then given by the formula (||). 

Now let us see that the extremum of W at the fixed energy e and n is in fact its maximum. Indeed, 
let us consider the functional 

' d d kb k b* k e uT 
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This functional is bounded below, so it has a minimum at some value of b], (at this point W as a function 
of bk has its maximum) , 

/ d d kb k b* k e uT = C{T) > const 

J min 

If one fixes some particular surface of singularities E, than the obtained value of the functional 



d d kb k b* k e ujT 



= C S (T) >C(T) 



for all values of T. Taking a saddle point value of 9, one gets 9 = — Inn + In C(T) and 

W(T) = n In n - n + ET - n In C{T) has an extremum at T\ , 

W(T)y: = n In n — n + ET — n In Cs (T) has an extremum at T 2 
Comparing W(Ti) and W(T 2 ) the following chains can be obtained 

W(Ti) > W(T 2 ) > W S (T 2 ) if Ti provides a maximum of W 

Wt,{T 2 ) < Ws(7i) <W(Tx) if T 2 provides a minimum of Ws 

Therefore, if one fixes some class of surfaces of singularities and energy e or, in other words, finds some 
particular singular solution, than one can get a lower bound for W(E,n). 

Due to the nonlinearity of the field equations it is not possible to obtain an analytical solution to 
the boundary value problem formulated above. One can find, however, analytical solutions in several 
regimes. 



3 Low energy expansion of F(An, e) 

In this section we calculate F(Xn,e) to the accuracy of 0(e 3 ). The typical momentum of final particles 
at small e is much smaller than the mass, so one can expect that the solution to the field equation ip c is 
a slowly varying function of x. The x-independent solution having singularity at r = and decaying at 
t — ► 00 can be found exactly Q, 

<Pq(t) = 

It corresponds, however, to the case of all final particles being at rest. Let us modify this solution by 
imposing by hand dependence on x. One writes 

^o(r,x) = y| sinh(r ^ ro(x)) 

This function has the surface of singularities r = to(x). In the end we will extremize the functional W, 
Eq.(fjj), with respect to ro(x). According to the general formalism described above, ro(0) must vanish, 
and due to spherical symmetry ro(x) depends only on |x|. Starting from this function (which satisfies 
the field equation to the accuracy of 0((9 x to) 2 )) one can expand cp in the following way, 

tp = ip + ipi + tp 2 + ■ ■ ■ (9) 

where <p\, ip 2 are of order of (9 x To) 2 , (<9 x to) 4 respectively. The explicit form of ip%, ip 2 can be found in 
Appendix. 

In order to obtain b k corresponding to the solution, let us find the asymptotics of (p at r — > 00 (see 
Eq.© l z = r-7D(x)) J 



1 

A sinh t 
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, , ,' z 2 4z 3 + tt 2 \ „ n/n n2 / z 2 z 5 \ ,„ 2 , 2 /z 2 2z tt 2 -7 n 
-9r 5 3 r — - — + — — + dT d{dT ) 2 __-+_+ (d 2 r ) 2 



12 J v ' I 2 2 12 / v ' V 4 3 12 

+ <^s(£-£ + ±) + ,^ 

This expression is of course the solution to the accuracy of O((9ro) 6 ) to the linearized field equation 
which has the following general solution, 



d d k bt 



Let us introduce the notation 

/" bl 



vwij&frjb*' (12) 



Thus, one gets, for any function of the form (fLl]), the following relation, 



r 2 



') - / ^s^-c - i* 2 + ^ **> - v> + ^ + <i3) 

Assuming that 

b = 2e T ° +6i+6 2 (14) 



where 61 ~ (9to) , 62 ~ (9tq) , one can find from ( |13[ ) that to the accuracy of 0((9tq) ), 

V = \/f e_T(2eT ° + ra2eT ° + 6l + ^~ 94eT ° + l d2bl + ^ ( 15 ) 
Comparing the terms of the same order in dr Q in Eqs.(|l5|) and ( px| ) one finds for 61 (0) and 62(0), 

h(0) = - 5 -9 2 ro, h(0) = - ^(^ro) 2 + ^Vr ) 2 (16) 



To get these expressions we took into account that ro(0) = and due to the 0(3) spherical symmetry, 
d z T (x)\ x = = 0. From Eqs.© and Q one finds the following system of the equations for b and its 
derivatives at x = 0, 

6 = 2- \d 2 r Q iaV„ - ^(a 2 t ) 2 + ^Vro) 2 ; 



a 2 6 = 29 2 ro - °-dl r d 2 T Q - -a 2 9 2 r , d% = 2d 2 r + 4(9 2 r ) 2 + 2(5 2 r ) 2 (17) 



11 

From this system one can easily find all the derivatives of tq with respect to x at x = via the derivatives 
of 6, 



d 2 T Q = \d 2 b+^b-^{d 2 3 b) 2 + \{d 2 b) 2 , (9 2 r ) 2 = i(9 2 6) 2 , 9V„ = i9 4 6-i(a 2 6) 2 -i(9 2 6) 2 (18) 
Substituting these solutions into the first equation of ( |l7j ) one gets the constraint on 6, 

2 - * + ^ + |> + ^44^ W + ^6^ ^ (19) 



From ( |12| ) and (19) one obtains 

\3ir 2 -25{[ d d k kikjb*^ 2 /A21-2tt 2 / f d d k k 2 b* k " ' 



4 144 \J (2n) d / 2 s/uj J V 4 96 \J (2n) d / 2 ^2 
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Let us now to extremize (Q) taking into account the constraint ( |2(i| ) . Making use of the Lagrange multiplier 
method, we rewrite (]?]) in the following form 

W = ET-n6- J d d kb k b* k e uT - e + A (^(b* k ) - y^J (21) 

where the constant A should be determined. Now W can be variated with respect to b\\ without any 
restrictions. Performing this procedure one gets 

A ^ T+6 (, 5 ,2 , 17,4 , , 21-2vr 2 



where 



C «=/(2^W> J <23) 
Assuming 6& = from ( p2| ) and fl23| ) one obtains the equation for Cy, 

r A fc^ T+e / 5 2 17 4 /A / 3vr 2 - 25 21 - 27T 2 2 

Uv ~ A J {2nY u> [ 12 fc+ 72 fc+ V4V 144 /c " fcm °" m + 96 ™ m 

(24) 

Let us now concentrate on the equation (24). First of all we consider the case i ^ j. Due to the fixed 
parity of the integrand, 

f - A A 3tt 2 -25 r d d k 2 2 T+g 
Gy -V4(27r)^ 36 Gy J {k]) e 

(no summations over repeated indices). From the last equation one concludes that CVj = 5ijTTC nm /d = 
5ijC/d. Together with Eq.(|3), this leads to 

C =(2^ e - T (l + °(i») (25 > 

T is assumed to be large at small e since, in fact, T ~ W eX tr/E. The constant A can then be found from 
the constraint ( |20| ) by making use of Eqs.(j2^) and (p5f): 

(26) 



fl6 e T ~ e / d(3d-26) d /9d 3 - 156d 2 + 884d - 1632 
A= VT(^r7 5 l 1_ 24T + 7^i x 1152 3tM 

where 

3tt 2 - 25 21 - 2tt 2 
72d + 48 

So, one finally obtains the expression for bk, 

IJq e T-uT T d/2 / d(3d-26) d /9d 3 -156d 2 + 884d-1632 
w V 24T + T 2 \ 1152 

The variation of the functional (pll) with respect to yields to 



_ A 2 e e ~ T / rf(3d-26) d_ f {d + 2)(3d 2 -58rf+ 272) \ 
" (27tT)^2 1/ + 24T + T 2 \ 384 + J 



(28) 
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From ( p8| ) 9 can be found as a function of T and n. The variation of (^l|) with respect to T gives the 
equation for T, 

d / 3d -26 13d - 102 8dA „ d 3d -26 / r 9d 2 -52d-140\ ,„„. 
^2t( 1 + ^2T IST^+T^ ^ T ^ + ^2- +£ ( 16 ' 72d ) ™ 

T is seen to be indeed large at small e. 

Collecting (Eq), (p8[), and (p9|) one finally gets for the function 



= n (in ^ - 1 + d - (in |- + l) + £ ^^1 JjLqtf 5 56d + 260 + 4 87 r 2 (d - 1))) (30) 

The first three terms in this expression coincide with those obtained in [0 and the last term is an e 2 
correction to W. Particularly, at d = 3, W has the form 

W** = n (in § - 1 + \ (in £ + l) - ,g + ^ (1327 - 9S. 2 )) (31) 

Finally, the last thing which is interesting to be performed is the calculation of To (a;) which corresponds 
to the obtained bk (gjj). To do it, one can find all derivatives of t (x) at x — (see (|l8|)) and after that 
restore To (a;) as the Taylor's series in powers of |x|. 

So, using the definition of b (|l2|), expression for (|27|), and taking into account that due to spherical 
symmetry dfj = %-9 2 one finds 

^ u +0 (^), ^.2^ +0 (^) (32 , 

All derivatives of To can be easily found from Eqs. ( [lS| ) and ({j2|), 

9 2 Tol „ = — — H -=77, 9? -Tn| „ = Sa | — — H — ) , 9 4 To| „ = dfii, m Tn\ _ = 9 2 A 2 -To| „ = 



Thus, the singularity surface which corresponds to the cross section up to 0(e 2 ) is a paraboloid, 

, / 1 5 1 \ x 2 £ / 3d - 16 \ 
^ = X (" 2T + 12 T 2 ") = [ l - £ ^ J < 33 > 

The Eq.(^3|) is valid for values of x being smaller than the typical inverse momentum of the final particles, 
e -1 / 2 . It is worth noting that the term of order (x 4 e 2 ) vanishes. These results are useful to compare with 
numerical calculations (see Sect. 5). 



4 Lower bound for F(Xn, e) in ultrarelativistic regime 

Let us now turn to the opposite limit, e > 1. In this case one can neglect the mass term in the field 
equation and consider the massless theory. Since Euclidean massless 4 theory is conformally invariant, 
several solutions can be found analytically. The simplest one is 0(4) symmetric Lipatov's instanton flifl . 
The corresponding solution which is singular at r = x = can be easily constructed Q by making the 
size of the instanton pure imaginary and shifting its center: 

P 11A\ 

V A x^ + (t + p) z - p l 

The variational procedure on T, 6 and only collective coordinate p (clearly, no collective coordinates 
describing the shape of the singularity surface appears) could be made and results in (see S) 

An 2 

W = nln An + n In — (35) 

16 7H 
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However, the variation on the shape of singularity surface is still required, unless the result (|35|) provides 
only lower bound for the true exponent. 

To see that these 0(4) symmetric calculations in fact do not provide the extremum over the shape of 
singularities, let us consider the small perturbations about the solution (p4|). Since the solution should be 
0(3) symmetric with respect to x rotations, it is sufficient to consider the approximate solutions which 
slightly vary the singularity surface so that the dependence of if), an angle between the r axis and the 
radius vector, appears. The solution should have the form 

Lp= — — — =cp + ^a m C^> (cos if>)<p m (r) + 0(a ) 

r 2 - p 2 {l + aG'in' (cos ip)) 2 m 

In this equation, (cosif)) are Gegenbauer polynomials, a m are small parameters, and the linear 
fluctuations (p m should satisfy the following equation, 

3 to(to + 2) 2 V 



To solve this equation, one can exploit the conformal invariance of the model and use the stereographic 
projection onto the sphere S (it is a convenient way to solve the equations in the conformal models, we 
use here a slightly modified method of Skipping the details, the solution is 

1 r 2 - m p m+1 ( 2(m - 1) p 2 (to - l)m p 4 

Vrn(r) = ^=—2 1 



y/X (r 2 - p 2 ) 2 \ to + 2 r 2 (to + 2) (to + 3) r 4 _ 
Calculating the asymptotics of this solution, one obtains the following expression for b^: 

~ - k of nm +ik m Q 2. (m+1) 



b k = bl + Y,*rJjTe- kp (p m+1 ^-3p 2 k o 
* — ' y U \ to! (m + 2)(TO + 3) 

(here 6° corresponds to the solution |||)). Then the exponent after the variational procedure on T and 
9 is modified by 



p V 2(to + 2)(to + 3) 

For to = or 1 (which corresponds to the dilatations and translations of (fo), AW = 0(a 2 ). However, 
for all other values of to the perturbation in the exponent is proportional to a m , so that the functional 
derivative of the exponent calculated at 0(4)-symmetric solution with respect to collective coordinates 
describing the shape of the singularity surface is non-zero. So, the result ( |35|) is indeed only lower bound 
on the tree cross sections in high-energy limit. One could expect that for lower energies (massive theory) 
this property will be conserved. 



5 Lower bound of F(An, e) at arbitrary e 

So, to obtain a lower bound for <jt TCC (E,n) at given E and n one could consider 0(4) symmetric field 
configurations |l6| ]. The singularity surfaces are then three-spheres in Euclidean space, and the only 
variational parameter is the radius of a sphere, p. Since the singularity surface should touch the origin, 
the sphere should be centered at 

= (~p,0) 

i.e., the 0(4) symmetric solution has the general form 



ip = <p(r), r = yjs. 2 + (r + p) 2 
At large r, the solution tends to the exponentially falling solution to the free field equation, 

Kx(x) e~ r 
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i.e., at r — > oo one has 



exp(-^x 2 + (r + p) 2 ) 



where the coefficient function A(p) is to be determined by solving the field equations under the condition 
that it has a singularity at r — p. From this asymptotics one finds 

b k = 2A{p)- 



'2u 

The "action" (Q) is then expressed through A(p), so W reads 

W = ET — nd — 8nA 2 (p)e- e Kl ^ P ~J ] (36) 

2p — T 

Now we extremize Eq.(^6|) with respect to p, T and 0. This leads to the following equations which 
determine the saddle point values of these three parameters, 

E = Aftf = K 2 (2p-T) 

n A(p) K^p-T) 1 ' 

l- f^y } (38) 

So, we look for the classical solutions which are singular at the spheres r 2 = p 2 and from their 
asymptotics obtain A(p), then express saddle point values of p, T and 9 through E and n (by making 
use of the Eqs. (p7|), (|38|)) and finally obtain the estimate for the exponent for the tree cross section (see 
Eq.(|)) 

W trcc = ET-n0-n 

It is straightforward to perform this calculation numerically for all e. The exponent for the cross 
section has the form 

F(Xn, e) = \n In ^ - An + Xnf(e) + 0(X 2 n 2 ) 
16 
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with the function /(e) plotted in Fig.l. At low energies our result matches the perturbative results, 
Eq.(|3l|). The fact that our variational approach leads to the exact results for Wtree at small e can be 
understood as follows. At small e, the curvature of the singularity surface is always large, and only this 
curvature is relevant for the evaluation of /(e) 0. In other words, the surface of singularities has the 
form (see Eq.(|3^)), 

r(x) = ax 2 + 0(x 4 ) + ■ • ■ 

and only the leading term is important at small e. Clearly, this leading term can be reproduced exactly 
in our 0(4) symmetric ansatz, and our result is exact at small e. At very high energies (small values of 
p) our field configuration tends to the 0(4) symmetric solution of the massless equations, Eq.(|4|). So, at 
e 1 it approaches the lower bound derived from this solution. The alternative lower bound on /(e) can 
be easily read out from ref. and it is also shown in Fig.l. This bound has been obtained by direct 
analysis of diagrams. As one can see from Fig.l, our bound is stronger than that of ref. Q. 



6 Conclusions 

To summarize, the tree multiparticle cross sections can be calculated using a semiclassical method based 
on singular solutions to classical field equations. At low kinetic energies e of particles in the final state the 
corresponding solutions can be obtained analytically and the leading exponent in the cross section can be 
found at least up to 0(e 2 ). Instead of proceeding to higher orders in e, one can solve the field equations 
numerically for different surfaces of singularity and after that apply the variational procedure at arbitrary 
e. For the simplest form of the singularity, the spherically symmetric one, the variational procedure over 
collective coordinates is incomplete and provides a lower bound on the tree cross section (the correct 
result at tree level should correspond to the true maximum of the exponent). To improve this estimate, 
one should consider more general field configurations. Beyond the tree level, the semiclassical Landau- 
like methods are still applicable, but the variational procedure is more complicated and corresponding 
solutions are known only in the simplest case of zero energies. 

The authors are indebted to V.Rubakov for numerous helpful discussions. We would like to thank 
M.Polyakov, P.Tinyakov and M.Voloshin for valuable discussions. This work is supported in part by ISF 
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The work of F.B. is supported in part by Soros fellowship for students, and the work of M.L. and S.T. is 
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Appendix 

The first two terms in the expansion (^|) have been obtained in |t| , 

2 f 1 1 /„, . cosh z — sinh z 1 coshz . .„ x9 coshz 



+ ^ V A i sinh^ + 2 l 9 T ° ( 3 ItoTz + Z ^~J + {d ^ Z ^r z ) j (39) 

(z = t — To(x)). The function ipi satisfies the following equation, 

d^ipi - (pi- SXiplvi = -df<Po = -dlipo{diTo) 2 + d T (p d?To (40) 
The solution ( |39| ) can be found with the help of the Green function, 

f = <A) + / daK(z,a)(d„<p d 2 T - a^ (9r ) 2 ) + C^f^z) + C 2 {x)f 2 {z) 



where 



,2cosh(z) /2"cosh(z) /sinh 2 (z) 3(r - r ) 
A(r) = -d^o = \l~ x -^y Mr) - yj \— A — + tanh(z) 



are the solutions of the homogeneous Eq.([l0]) and 



K{z,a) = ^(f 2 {z)h{a)-f 2 {a)h{z)) 
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is the Green function. The functions C\(x) and C^{x) can be found from the conditions (p(r — > oo, x) — > 
and ip{0, 0) — > oo. 

The next correction (we expect that it is of order of (c? To)) can be obtained in the same way from 
the following equation, 

d?ip 2 - ¥>2 - 3A< y 2 ; </?2 = -d?(pi + 3\tpl(p 

and has the following form, 



[2 ( _o n2 / 1 1 1 , , , x 1 , . . . 1 z 1 cosh z 

- 1 t { — « y.Tn : (coshz — sinnz zlcoshz — smhz A z ^ — 

^ VA\ 1 3 V 8smh2 8 12 4 sinh z 8 sinh 2 z 

l 2 coshz\ & (a ,2(1. 1 1 cosh* 1 , _ , _,A , „ / 1 1 



— o- - dudiT y —— — -z — 5 -zcoshz + i(z) + 2aa 2 T • a ? T — -r-r— 

8 sinh 2 zy jV ; V 12 sinhz 3 sinh 2 z 12 V 7 Jj jU Vl2sinhz 

1 coshz 1 z 1 9 coshz 1 . , ,, A „„ r ,„ „ /l 1 1 coshz 

-z 2 - — 3 — — zsinhz + /(z) +2a J [(a i T ) 2 ] -a-T --r-r — o^" 



3 sinh 2 z 4 sinh z 4 sinh 2 z 12 / J J \8sinhz 8 sinh 2 z 

l/i • i \ 1 2 cosn z \^2^2/l , 1 • i l/i • i a 3 1 

- — (cosh z — smh zj + — z — — „- — + "i T o^i T o ~zz cosh z — — z sinh z + —(cosh z — smh z) — ——r~, — 
24 8 sinh z J \12 4 4 8 smh z 

3 coshz 1 z 1 z 2 1 z 2 1 2 coshz A 2 2 /l . , 

-z o ; 1 5 1 ; 1 — z o h27(z + (OiTn) o^Tn — z ( cosh z — sinh z) 

8 sinh 2 z 2sinhz 4 sinh 3 z 8sinhz 4 sinh 2 z KJ J V ' 3 \12 v ' 

1 , , . , , 1 z 2 1 z 2 11 3 coshz 



-(cosh z — sinh z) H = 1 1 z ., 

12 2 sinh 3 z 4 sinh z 4 sinh z 4 sinh 2 



where 



1 i /r. • i \ / • i 3 coshz \ 1 coshz , 

i(z) = — In 2smhz) smhzH 3z o — A ^ — / xcothxefcc. 

12 v A sinhz sinh 2 z 7 4 sinh 2 ' 
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